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There have been many proposed methods for the prac- 
tical implementation of quantum computing. Now quantum 
computation has reached the turning point from being a con- 
ceptual system to becoming a physical one. In this paper, we 
discuss a practical elementary gate and the actual computa- 
, tional time-cost of the QFT in two physical implementations, 
' namely the bulk spin resonance computer and the Spin Res- 
onance Transistor. We show that almost all universal gates 
require different times for operation. The actual time-cost of 
^ the QFT is 0(n2") for large n. This differs drastically from 
, the reported cost O(n^) based on ideal quantum computation. 



quantum computers being considered are constructed 
from an array of n quantum registers labeled by j(l < 
i < n) in one dimension. They calculate n qubit data 
which corresponds to = 2" states. We defined the fol- 
lowing matrices. The suffixes on each matrix represent 
the quantum registers on which it operates. In the def- 
initions of Cj^k{d) and Dj_k.{d), the first and the second 
suffixes represent the target bit and the controlled bit 
respectively. 
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Recent technological development has stimulated pro- 
posals for many quantum computers: Bulk Spin Res- 
onance(BSR) (f]-|], trapped ions g, cavityQED 
Josephson junctions ||^, coupled quantum dots and 
the Spin Resonance Transistor(SRT) BSR has been 

implemented experimentally in some organic molecules 
by the use of conventional nuclear magnetic reso- 
nance (NMR) equipment the SRT is attractive from 
the viewpoint of it's integration ability and it's compat- 
ibility with silicon technology. 

Complexity analysis classifies quantum algorithms ac- 
cording to a function that describes how a computational 
cost incurred in solving a problem scales up as larger 
problems are considered The computational cost of 
a quantum algorithm has usually been estimated as the 
sum of the universal gates required in such ideal mathe- 
matical models as the Quantum Turing Machine(QTM) 
and the quantum circuit. The computational complexity 
is effective in estimating the essential performance of an 
algorithm to factor out the variations in performance ex- 
perienced by different makes of computers with different 
amounts of Random Access Memory (RAM), swap space, 
and processor speeds. The above cost is proportional to 
an actual time-cost in the physical implementation where 
all quantum operations can be achieved in the same time. 
However, if the implementation being considered takes a 
different time for each quantum gate, there is a possibil- 
ity that the actual time-cost will have a different behavior 
from the ideal cost. A hardware dependent time-cost is 
important to experimentalists who research into practical 
implementations. In this paper, we focus on the actual 
computational time-cost of the Quantum Fourier Trans- 
form (QFT) in two physical implementations that utilise 
nuclear spins: BSR and the SRT. 

For the following discussion, we assumed that the 
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Firstly, we discuss the practical elementary gate in 
these implementations. The practical elementary gate is 
defined as a quantum gate that can be achieved directly 
by physical phenomena in the implementation being con- 
sidered. 

BSR consists of quantum registers that are the nuclear 
spins of each atom in an organic molecule . External 
RF-pulses and magnetic fields control the quantum states 
and the quantum correlations in the following way. 

The nuclear spin under a strong magnetic field 
B = (0, 0, So) is described by the Hamiltonian R = 
—"fhBoIjz, where 7 is the gyro- magnetic ratio for the 
spin, and Ijz is the z component of the j'th nuclear 
spin. The time evolution of the system is exp(iHt/h) = 
exp{—i"/hBotSjz) and this then corresponds to the rota- 
tion on the z-axis Rzj {0 — jhBotY The RF-pulse enables 
the rotation on the other axis ||ll|. The time evolution 
shows that the rotation angle 6 can be controlled by two 
factors: the intensity -Bo, and the duration t, of the exter- 
nal magnetic field. Consequently, there are two control 
modes. In the intensity control mode, each single qubit 
rotation takes the same time. In the duration control 
mode however, each single qubit phase rotation takes a 
different time, which is proportional to the rotation an- 
gle. 

The exchange interaction between the j-th and the k- 
th registers is described by the Hamiltonian Hf-^ch. = 



1 



Jjkljzlkz, where Jjk is the exchange coupHng constant. 
The time evolution of the system is described by 



exp{iJjktIjzIkz/K) = Dj^k{d = Jjkt/2h). 



(1) 



Sequence(|^) shows that the exchange interaction can 
control the phase rotation angle. Any external fields, 
however, cannot control the interaction directly, because 
atoms in a molecule always interact with each other. The 
refocusing technique can control the angle effectively 
[ pl| . Only the time duration between the refocusing 
pulses determines the angle 9. The operation Dj k{0), 
therefore, takes a time that is proportional to the rota- 
tion angle 9. An actual operation such as Dj,k is effective 
only for adjacent registers because the exchange interac- 
tion Jj^k between non-adjacent registers is very small. 

The SRT is composed of quantum registers which are 
the nuclear spins of arrayed phosphorus ions in silicon, 
with a globally static magnetic field B and an AC mag- 
netic field Bac § • The implementation consists of two 
gates on the surface: the A-gate above each ion and the 
J-gate between adjacent ions. 

The A-gate controls the strength of hyperfine inter- 
actions and the resonance frequency of the nuclear spin 
beneath it. A globally applied magnetic field Bac Aips 
nuclear spins resonant with the field by the same process 
that occurs in BSR. In this case, only the duration of the 
resonance determines the rotation angle. It implies that 
each single qubit phase rotation always takes a different 
time. 

The electron wave function extends over a large dis- 
tance and makes an effective electron-mediated coupling 
for two nuclear spins sharing it in semiconductors. The 
J-gate controls the overlap of the electron wave functions 
bounded to two adjacent phosphorus atoms, and hence 
the electron-mediated exchange coupling Jj^k — J{t) in 
the time evolution(|l|) directly. It therefore controls the 
phase rotation angle in the operation Dj ^- The opera- 
tion Dj^k in this implementation also operates only for 
adjacent registers. 

In both implementations, all single qubit phase rota- 
tions and controlled phase rotations are practical elemen- 
tary gates. They can make the quantum XOR in the 
sequence below, ordered from right to left 
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The suffixes j and k represent the target bit and the 
controlled bit respectively. The sequence (||) shows that 
the quantum XOR depends on single qubit gates in this 
technique, and is then concerned with the time required 
for phase rotation. 

Barenco and co-workers showed that other universal 
gates can be constructed by all single qubit gates and 
the quantum XOR iQ. Each n{> 1) qubit gate required 
a different number of them fl^ , and then a different time 



for execution. As discussed above, almost all univer- 
sal gates take a different time in these implementations. 
It suggests the possibility that the actual time-cost of 
the quantum algorithm is different from the ideal cost, 
though the complexity is not affected in the QFT case. 

Next we estimate the actual time-cost of the QFT 
that could be achieved by the above practical elementary 
gates, by considering the time resolutions of the control- 
ling external fields in these implementations. 

The QFT is the transform with base N — 2", corre- 
sponding to the n qubit defined by 
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Shor proposed the algorithm factoring a composite inte- 
ger by the QFT H. The QFT can be constructed by 
the sequence in the order (from right to left) 

-f^oCo,l(^l)C'o, 2(^*2) • • • Co^n~l{6n-l)Hi ■ ■ ■ i/„_3 
Cn~?,,n-2{0l)Cn~3,n-l{92)Hn-2Cn-2,n-l{9l)Hn-l (4) 

followed by a bit reversal transformation, where 9j = 
. The n qubit QFT requires n{n — l)/2 controlled 
phase-shifter {C'j^ki9k~j)s and n Hadamard transforma- 
tion HjS, and then n{n +l)/2 ~ 0{n^) This esti- 
mation is based on the complexity analysis method. It 
coincides with an actual time-cost in the case where all 
gates required in the sequence(^) are practical elemen- 
tary ones. 

The controlled phase rotation Cj^k{9k-j) can be 
achieved by the sequence 

Cj,k{9k-j) ~ Rzk{—9k-j+i)^k{9k-j+2)Rzj{—9k-j+i) 
XORU, k)R,j{9k^j+i)X0R{j, k). (5) 

The intensity control mode can make all operations 
{Cj,k{9k-j), Hj) take almost the same time. Conse- 
quently, the actual time-cost coincides with the ideal cost 
in this mode. 

The duration control mode, however, makes each 
operation Cj^k{9k-j) require the time Tk^j in propor- 
tion to the phase rotation angle 9k- j. The operation 
Cf)^n-i{9n-i) rotates the minimum phase angle 9n^i and 
takes the minimum time t„_i of all rotations in the 
QFT(4). The range of required phase rotations in the 
QFT(4) increases with 2". For example, the time ratio 
tq/tVi-i approximates to 2-'^'"' ~ 10^° in the 100 qubit 
QFT. 

In general, the time resolution controlling the exter- 
nal field is determined by the response time of the system, 
the delay of the electronic signal and so on. We can only 
execute in the physical implementations that satisfy the 
relationship: 



To > Tl > • • • > T„_l > tji. 



(6) 



It is important for the actual time-cost estimation to 
determine how we set up the unit time tumt- The unit 
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time tunit should be also greater than the time resolution 
t/j. The QFT's in these implementations have various 
actual time-costs from tunit — tq to tunit — Tn-i- 

If we adopt the maximum rotation time tq as the unit 
time tunit, the actual time-cost is 0{n), since 

n—2 n— 1 n—2 n— 1 ^ n—2 n— 1 

Tk-j _ \ ^ \ ' tik-j _ \ ^ \ ^ 
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= n + 2i-"-2~0(n). (7) 



On the other hand, the condition tunit 
the actual time-cost 0(n2"), since 



= r„_i makes 
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j=0 k=j + l 

= (n-2)2"-i + l~0(n2"). 



(8) 



In this way, the actual time-cost varies from 0(n) to 
0(n2"), and it depends on which of these is adopted as 
the unit time tunit- The former time-cost however, is not 
valid for any n in the following way. 

The former condition, tunit = tq, means that all phases 
are always rotated by the external field with constant in- 
tensity B for various data of magnitude n. In this case, 
the minimum time t„_i decreases exponentially with in- 
creasing n. We cannot rotate the phase 6j to satisfy the 
condition Tj < t^. There exists the upper bound ni,, 
satisfying the relationship(^ for the intensity B under 
consideration. The estimated time-cost(|^) is valid for 
any n satisfying n < rif,. It is, however, not valid for any 
n which is greater than rtf,. 

The latter condition tunit = Tn-i means that the in- 
tensity B decreases exponentially with n satisfying the 
relation(|^). In this case, we can rotate all phase angles 
^i(0 < j < 1) in the QFT(^ accurately. A particular 
condition tunit = Tn-i — til always achieves all phase 
rotations in the minimum total time, and then yields the 
best computing performance for each value of n. The 
actual time-cost always obeys eq.(^ for any n. 

In this way, the actual time-cost varies from 0(n) to 
0(n2") for any n(< rif,), and follows only the latter for all 
other values of n. These costs are estimated making the 
assumption that the QFT is always executed accurately. 

An approximate QFT, (AQFT) can reduce the arbi- 
trary numbers of the controlled phase shift gates by sac- 
rificing the accuracy We can select the AQFT with 
an actual time-cost between 0{n) and 0(n2"), by con- 
sidering the required accuracy for any n. 

We have discussed the actual time-cost from the view- 
point of the phase rotations in the QFT(^. Almost 
Cj^k{Qk-j) operations occur for non-adjacent registers in 
the QFT. We need to construct such non-adjacent gates 



using adjacent ones so here we estimate the actual time- 
cost required for constructing such non- adjacent gates 
from adjacent ones. We must construct any non-adjacent 
gates by use of the adjacent swap technique. The swap 
Sj^k is an operation for exchanging data between two 
quantum registers simply. It is achieved via the se- 
quence 5j-fc = XORik,j)XOR{j,k)XOR{k,j). The 
non-adjacent two qubit gate Uj,k is achieved by adjacent 
swaps and adjacent C/;.i+i through the following fc — j — 1 
sequences. 

Uj^k — Sk^k-lSk-l,k~2 • • ■ Sj+z,j+2Sj+2,j + lUjJ+l 

^j+l,j+2Sj+2,j+3 • • • Sk~2,k-lSk-l,k (9) 
= Sk,k~l ■ ■ ■ Sl+2,l + lSjJ+l ■ ■ ■ Sl-ijUl,l + i 

■ ■ ■ 'S'j+ij-S'/+i,;+2 • • • Sk-i,k{j < I < k) 

^ ^j-j + i^j + i-j+2^j+2-j+3 ■ ■ ■ Sk-2,k-lUk-l,k 

Sk-l,k-2 ■ ■ ■ Sj+3j+2Sj+2,j + lSj + l.j. (10) 

Figure 1 represents the sequence^. These sequences 
show that the non-adjacent operation Uj^k requires both 
the time for the adjacent operation Ujj^i and another 
time in proportion to |fc — j| for swaps transferring data. 

Considering the construction dgj^-^), the QFT(^ re- 
quires more {n — l)n{2n — l)/6 ~ 0{n^) adjacent swaps 
for data transfer besides the adjacent controlled phase 
rotations. 

There is, however, a way to reduce some swaps. 
From sequences (^) and (|l0|), we can obtain the actual 
quantum circuit shown by Fig. 2. In the figure, swaps 



^3,j+i ' ■ ■ Sk-l,kSk,k-l ■ ■ ■ 



inside each box can be 



reduced to the identity operation/. The n qubit QFT 
requires {n — l){n — 2)^ 0{'n?) swaps as a result. It has 
the same polynomial order as the ideal cost. In this way, 
there are some cases where the swaps can be reduced due 
to the symmetry of the algorithm under consideration. 

We can conclude that the actual time-cost of the QFT 
is dominated by the phase rotations and then are 0(n2") 
in the range n > rih. 

The range of required phase rotation angles increases 
exponentially with n in the QFT if the accuracy is pre- 
served. In the implementations that are considering, we 
need to obtain it by controlling the duration or the in- 
tensity of the external field. The duration control mode 
increases the actual time-cost drastically with 2". The in- 
tensity control mode in BSR takes least time and seems 
to be the most efficient case for the QFT. This mode, 
however, leads to another burden on the equipment. The 
required range of intensity for the applied field increases 
with 2". If the minimum phase rotation is implemented 
by the field B = 10~^ T, the maximum phase rotation 
for the 100 qubit QFT requires a field intensity B 10^^ 
T, which is far beyond the current feasible intensity for 
the magnetic field. 

Our results lead to concerns about the feasibility of 
factoring huge numbers with polynomial order time-costs 
in the implementations we are considering. 
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We have shown that almost all universal gates do not 
expend the same time in BSR and the SRT. This causes 
the actual time-cost of the QFT to be drastically different 
from the ideal one. The ideal cost of the QFT is not 
effective in the range n > rib for practical cases of these 
implementations if the accuracy is preserved. 

We believe that both ideal and actual discussions are 
important for the development of quantum computer sci- 
ence. The former discussion stimulates study of the char- 
acteristic of the algorithm itself, and the latter is also 
important for the current situation in quantum compu- 
tation as we move from the conceptual system to the 
physical one. Our discussion shows the necessity of dis- 
cussing the practical elementary gate in other proposed 
quantum computers, and of estimating the actual time- 
cost for other quantum algorithms. 
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FIG. 1. Non-adjacent operation achieved by utilising ad- 
jacent operations. The symbol inside the box represents the 
adjacent swap Sj^k- 
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FIG. 2. Figures (a),(b) and (c) in order are the actual quan- 
tum circuits of the 5 qubit QFT. This quantum circuit is de- 
rived from sequences (4) and (10). The swaps inside the boxes 
are reducible. 
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